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Abstract 

In this series of papers we examine magnetic reconnection in a domain where the magnetic field does not 
vanish and the non-ideal region is localised in space. In a previous paper we presented a technique for 
obtaining analytical solutions to the stationary resistive MHD equations in such a situation and examined 
specific examples of non-ideal reconnective solutions. Here we further develop the model, noting that 
certain ideal solutions may be superimposed onto the fundamental non-ideal solutions and examining the 
effect of imposing various such flows. Significant implications are found for the evolution of magnetic fiux 
in the reconnection process. It is shown that, in contrast to the two-dimensional case, in three-dimensions 
there is a very wide variety of physically different steady reconnection solutions. 

1 Introduction 

Magnetic reconnection is a fundamental process in astrophysical, space and laboratory plasmas. The first 
models of the process were two-dimensional (2D) and stationary; much of our present understanding of the 
subject (for a review see, for example, Priest & Forbes, 2000) relies on these models. However, astrophysical 
plasmas systems are typically three-dimensional (3D) with any non-ideal regions localized in each dimension. 
Models incorporating these features have demonstrated important differences between the 3D and 2D cases 
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(Hornig & Priest, 2003). In this paper we focus on the basic properties of 3D steady magnetic reconnection 
in a resistive MHD plasma with the aim of exploring these new features. 

In 2D, a change in the magnetic connectivity of plasma elements occurs at an X-type magnetic null- 
point that is located in a non-ideal region of an otherwise ideal plasma. Through analytical modelling, a 
range of steady-state 2D resistive MHD reconnection theories have been developed (e.g., Petschek, 1964; 
Vasyliunas, 1975; Priest & Forbes, 1986; Priest & Lee, 1990; Craig & Henton, 1995). Expansion procedures 
are one of the techniques that have proved invaluable in the development of these models, with the small 
parameter commonly taken as the Alfven Mach number of the flow. The early paper by Petsheck (1964) 
used such an expansion and, subsequently. Priest & Forbes (1986) applied the method to develop a family of 
almost-uniform regimes in which the nature of the inflow velocity determining the particular regime and with 
Petschek's mechanism being a special case. However, all 2D models have a number of elements in common; 
a stagnation flow that crosses the separatriccs of the field brings magnetic fiux in toward and, subsequently, 
away from, the non- ideal region and, with the electric field being uniform, the rate of reconnection is governed 
by the solution in the ideal region. 

In 3D, reconnection may also occur at null points, now in various forms depending on the direction of 
the current (Pontin et al. 2004, 2005a), as well as along separators connecting null points (Lau & Finn, 
1990, Priest & Titov, 1996) and at locations where the field does not vanish (Schindler et al. 1988, Priest & 
Forbes, 1992). This non-null reconnection requires only a localised non-ideal region together with a non-zero 
integrated electric field component parallel to the magnetic field (Hesse & Schindler 1988). Each of these 
forms of 3D reconnection may look quite different, both from each other and from the heavily analysed 2D 
case. 

Regarding the specific case of non-null magnetic reconnection and a localised non-ideal region, the sub- 
ject of the present study, Hornig & Priest (2003) presented a kinematic stationary model of the process. 
Here the term kinematic refers to the situation where the magnetic field is imposed and the equation of 
motion neglected. The authors noticed that in 3D, under the kinematic assumption. Ohm's law can now be 
decomposed into non-ideal and ideal components: 

"^non— ideal ~t~ ^non— ideal ^ B = TfJ 
^ideal + ^ideal X B = 

As a consequence of the localisation of the non-ideal region, the plasma flows in the non-ideal solution, termed 
the 'particular solution' are rotational in nature, with the rotation being in different senses above and below 
the non-ideal region. An important feature of the particular solution is that it is spatially restricted with the 
rotational flows being found only within the field lines threading the non-ideal region. Onto this particular 
solution can be imposed any solution to the ideal Ohm's law, the results being termed 'composite solutions'. 
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Such solutions can ease the spatial confinement of the reconnection process, allowing its effect to be seen on 
a global scale. However, the independence of the non-ideal and ideal solutions suggests the characteristic 
coupling in 2D between large-scale flows and reconnection rate no longer holds in 3D. 

In a previous paper (Wilmot-Smith, Hornig & Priest, 2006, hereafter Paper I) we introduced an expansion 
scheme through which we may carry out a dynamic analysis of 3D non-null reconnection (i.e. including also 
the equation of motion) and so address fundamental questions regarding both the apparent freedom in the 
3D system (as suggested by Hornig & Priest, 2003), and the nature and diversity of non-null reconnection 
solutions. In that introductory paper we developed the solution method and demonstrated how the system 
can be broken down into particular and composite solutions. We then proceeded to examine the case 
of particular solutions, again finding counter-rotational plasma flows confined to be within the flux-tube 
threading the non-ideal region. However, in the dynamic analysis the ideal and non-ideal components to 
the 3D Ohm's law are coupled via the equation of motion. Although we may again impose an ideal solution 
onto the particular solution and create composite solutions, the coupling will restrict which ideal solutions 
may be taken. We address this important question in Section [3] and find certain restrictions on the choice 
of ideal flow. Accordingly we then proceed to examine an example of a more general ideal flow in Section |3] 
where the precise division between particular and composite solutions no longer holds but, nevertheless, a 
great deal of freedom within the 3D system is found. 

2 Model 

Using a method described in more detail in Paper I, we consider the dimensionless, stationary, incompressible, 
resistive MHD equations and seek three-dimensional solutions for slow flows, {v ^ va)- By assuming the 
inertial term in the momentum equation is small, we may expand the dimensionless variables in terms of 
the Alfven Mach number of the flow, <C 1: 

B = Bo -I- MeBi + M^B2 + M^Bs + • • • , 

V = Vi + MeV2 + M^v-s + ■■■ , 

j = jo + A/ejl + Me j2 + Me ja + ' ' ' , 

E = Eo + 7\/eEi + Af |E2 + • • • 

= -V>o - A/eVVi - M2V'02 + • • • , 
p = PO + MePl + MIp2 H , 

where all the quantities B^, v^, j;, E;, (pi, pi are dimensionless. 
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As in the analysis of Paper I we assume throughout that the lowest-order magnetic field, Bo is potential 
(and so jo = 0) and that Bi = 0. With these assumptions a comparison of the zeroth- and first-order terms 
in the expansion of Ohm's law yields ideal and non-ideal equations, respectively. These are given by 

Eo -I- vi X Bo = 0, (1) 

E1+V2 xBo = 77j2. (2) 

In order to satisfy the continuity equation we are assuming for simplicity that the plasma is incompressible. 
Accordingly, V • vi = V • V2 = 0, together with ([T]) and are the only non-trivial equations of zeroth and 
first order. The pressure terms po and pi are both constants, and the main dynamic effects are described by 
the equation of motion at second-order: 

(vi-V)vi=-Vp2+j2xBo. (3) 

This equation provides a coupling between ([T]) and 
Throughout the analysis we set 

Bo = bo{ky,kx,l), (4) 

where fc > 0, so that the basic state is a current-free equilibrium representing the superposition of an X-type 
null-point field in the xy-plane with a uniform z-component of the magnetic field. The equations, X(xo, s), 
of the field line passing through the point xo are given by 

X = cosh(&o/cs) -f yo sinh(6ofcs), Y = i/q coah{bQks) + xq sm]i{bQks), Z = bo-s + zq, (5) 

with the inverse mapping Xo(x, s) given by 

Xo = .T cosh(6ofcs) — J/ sinh(5oA:s), Yo = y cosh(&ofcs) — x sinh(6o^s), Zq = —bo-s + z, (6) 

where the parameter s is related to the distance I along a field line by ds — d//|_Bo|- 

The lowest-order equations in the expansion procedure allow for a direct comparison with the particular 
and composite solutions of Hornig & Priest (2003), whose 'particular solutions' are recovered in this scheme 
by setting Eo = and vi = 0. These solutions were described in some detail in paper I, where two examples 
of analytical solutions to ^ together with ^ were given. In both cases the non-ideal term, 1)^2, is localised 
in all three-dimensions. As a result of that localisation the reconnective plasma fiow, V2 is found to be 
confined to within the HPT and counter-rotational in nature i.e. rotating in opposite senses above and 
below the central xy-plane. Thus it is only the magnetic flux linking the non-ideal region that is affected by 
the reconnection process. 

In realistic situations we might expect a non-zero plasma velocity outside of the HPT. One way to achieve 
this is to take a non-zero ideal plasma flow vi. Whilst in the kinematic analysis an arbitrary ideal flow may 
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be superimposed on the particular solution, in this dynamic analysis the zeroth-order (ideal) and first-order 
(non-ideal) Ohm's laws HI [2]) are coupled through the equation of motion (in particular equation [3]). In the 
following section, Sec. [31 we describe the circumstances in which an ideal flow may be superimposed onto 
the same solutions for ^ that were detailed in paper I. We describe the nature of these ideal flows and the 
consequences of their inclusion on the evolution of magnetic flux. Certain restrictions on the choice of ideal 
flow that may be superimposed on the particular solution are found and so we in proceed, in Section 01 to 
consider a more general ideal flow with a slightly different physical motivation. 

3 Composite Reconnection Solutions 

In general, the momentum equation given by ^ implies a coupling between the ideal and non-ideal Ohm's 
laws given by equations (H]) and However, for the class of ideal plasma flows vi for which the curl of the 
inertial term on the left hand side of ^ vanishes, the equations become decoupled. In this case the effects 
of a non-trivial solution to ([1]) are apparent at second order only in the form of the pressure gradient, Vp2- 
For ideal flows satisfying this condition the particular solutions of paper I may be taken as a solution to ^ , 
and so we have a direct comparison with the composite solutions of Hornig & Priest (2003). 

We begin in Section 13.11 by examining an ideal stagnation flow vi for which V x (vi • V) Vi — and 
continue in Section [3. 21 to consider the implications on the evolution of magnetic flux on both this stagnation 
flow and other ideal flows with curl-free inertial terms. 

3.1 Simple Stagnation Flow 

Clearly there are a wide variety of ideal plasma flows vi for which V x (vi • V) vi =0 and are therefore 
candidates for the formation of composite solutions. The different ideal flows each correspond to physically 
distinct scenarios. Stagnation flows are an obvious choice to consider since they enable the formation of thin 
current sheets. In addition they allow for magnetic flux to be brought into and removed from the HFT and 
so, in contrast to the particular reconnection solutions, for the effects of the reconnection process to be seen 
on a wider scale. 

We now turn to equation [T] and determine a suitable ideal solution. We let 0o be the function of the 
field-line coordinates given by (xq , yo ) given by 

A ^0 !'7\ 

90 = j^xoyo- (7) 

Then, setting Zq = the inverse field line mappings © can be used to find an equivalent expression in terms 
of X, y, and z, so determining (f>o{x,y, z). The component of vi perpendicular to Bp may then be deduced 
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(a) (b) 

Figure 1: The ideal plasma velocity vi for (a) z = 0.5 and (b) z — —0.5, and the parameters (^q = 1, fc = 0.5, 
6o = 2, and A = 1 



from ([T]) as 

V^o_x Bo 
|Bo 



Vl-L = IP 12 ' 



and the freedom in choosing a component parallel to Bq exploited to set the z-component of vi to zero and 
so ensure the flow is divergence-free: 

(vi^)^Bo , - 

vi=vi^ . 9 

Thus we obtain an ideal stagnation-type flow given by 

Vi = (a;cosh(2fcz) — y sinh(2fcz)) ;e H — (xsinh(2fcz) — y cosh(2A;z)) y, 

for which V x (vi • V) vi = 0, as required. Some streamlines of the flow vi are illustrated for a particular 
choice of parameter values in Figure [TJ 

Since the inertial term (vi • V) vi may be expressed as the gradient of a scalar function, the equation of 
motion [3] has the same structure as that in the case of the particular solutions examined in Paper I (where 
vi =0). We therefore take the expressions for j2, ?), Ei and V2 that were given in Paper I, Section 3.2. The 
effect introducing the ideal-flow is seen in the pressure term p2 which now becomes 



V2 = P20 - |fcA 5oj20 tanh j - -j^ . (10) 

The corresponding term for the particular solution may be recovered by setting ifQ = 0, and so it is seen 
that the inclusion of a zeroth-order flow has had the effect of introducing an extra term to the pressure, 
proportional to (^q/ (A'^^q). When (po = Q there are strong gradients in the pressure along the separatrices of 
Bq in the xy plane. This extra term has the effect of smoothing out these strong gradients, with p2 becoming 
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Figure 2: The pressure profile p2 [x, y) when vi 7^ 0, j2 = ^20/ cosh^ ((x^ — j/^) /A^) and the parameters 
P20 — 2, bo — 2, X ^ 1, k = 1/2, ^20 = A = 1 and (po — 1/2. The lower pressure regions correspond to 
inflow of vi and the higher pressure regions to outflow. 



a smoother function as Lp\l l\^h\ is increased. An example of the pressure profile is shown in Figure^ which 
can be compared with Figure 6 of Paper I. 

The additional term has a natural physical explanation. It deflects the incoming vi flow toward the 
outflow direction, a purely hydrodynamical effect. Due to the symmetry with respect to inflow and outflow, 
there is no net transfer of magnetic energy to kinetic (bulk) energy of the plasma in this stationary solution, 
as would be expected in a more realistic situation. However, we may model part of this process by requiring 
vi • j2 X Bo to be positive. This would result in an initial transfer of magnetic energy to kinetic energy, but 
with the latter subsequently transferred to potential energy, since v • Vp > 0, so no net acceleration can take 
place. We have here that 



In order to model the situation in which the Lorentz forces drive the reconnecting plasma flow we require the 
expression (jlip to be positive, which can be ensured by taking the combination ipQj2nk < 0. Note, however, 
that cases where external pressure differences drive flows against Lorentz forces ('driven-reconnection') are 



vi • j2 X Bo 



Vahok {{x^ + y'^) cosh {2kz) — 2a;y sinh (2/cz)) 



(11) 
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also known. 

The reconnection rate in the model is given by the integral of the parallel electric field along the recon- 
nection line which is identified with the ^;-axis: 

dt y II ' J |Bp| 

= v^Me eipL + O (M2) . (12) 

Since we have taken the same non-ideal electric field Ei in both the particular solutions of Paper I and these 
composite solutions the expression for the rate of reconnection is clearly the same in both cases. However, the 
non- vanishing external ideal flow changes the interpretation of this reconnection rate since the reconnection 
process can now reconnect flux lying initially outside of the HFT, as will be demonstrated in the following 
section. 

We could, in principle, determine higher-order quantities of the expansion solution using the iterative 
scheme outlined in Paper I. In the particular solution wc found that Ohm's law at subsequent even-orders 
and the equation of motion at subsequent odd-orders could be set to zero but in this case the remaining 
equations must all be solved numerically. 

3.2 Magnetic Flux Evolution 

We have shown that the introduction of the ideal stagnation flow vi, as described in the previous section, 
does not alter the reconnection rate (up to second order). The evolution of magnetic flux in the two cases 
is, however, quite different, and may be visualised using the concept of a magnetic flux velocity (e.g. Hornig 
& Schindler, 1996). Such a velocity satisfies 

E-l- w X B = 0. 

In two-dimensions, a non-vanishing electric field at an X-point of the magnetic field requires a singularity 
in the flux velocity at that point. The singularity describes the reconnection taking place at that point. 
In 3D, on the other hand, there is no general well-defined flux-transport velocity. However, provided the 
non-ideal region is localised and contains no closed flux, we can replace the notion of a single flux-transport 
velocity by a pair of velocities Wi„ and Wgut that describe the behaviour of magnetic flux entering and 
leaving the non-ideal region, respectively (Priest, Hornig & Pontin, 2003). The two flux velocities will not 
coincide within the non-ideal region and through examining the difference between the two we may deduce 
information about the behaviour of magnetic flux in the process. 

The use of these flux velocities therefore allows us to track the evolution of plasma elements in the 
ideal region above and below the localised non-ideal region. Initially-connected elements will only remain 
connected if the field line linking them does not pass through the non-ideal region; otherwise the elements 
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Figure 3: Flow Wi„ (grey) and Wo„t (black) for the solution described in Section [3{a). 
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X 

Figure 4: Separatrices of Wi„ (grey) and Wout black for the solution described in Section [SJa). The region is 
divided into three types of reconnective behaviour. Magnetic flux passing through region I undergoes ideal 
evolution. Magnetic flux passing through region 11 undergoes a slippage-like behaviour while flux passing 
through region II undergoes an evolution similar to that seen in classical 2D reconnection 
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will change their connectivity. The flow-lines corresponding to the ideal evolution above and below the 
non-ideal region can be projected onto the central plane, say, using the pair of quasi-velocities Wi„ and 
Wout- For the stagnation flow described in the previous section, the relevant projection is shown in Figure [3] 
(for a particular choice of parameter values). The flow lines of Wm (grey lines) in the z — plane are 
superimposed on those of (black lines) in the same plane. We are able to divide the plane into three 
regions according to the type of reconnective behaviour that occurs; the separatrices dividing these regions 
are shown in Figured] 

In region I the flow lines of Wi„ and Wo^t coincide perfectly. The magnetic flux passing through the 
z — plane in region I evolves ideally, so that initially-connected plasma elements will remain connected. 
In regions II and III the flow lines of w,;„ and "Wout do not coincide. For magnetic flux passing through 
the 2 = plane in these two regions we deduce that plasma elements above and below the non-ideal region 
that are initially connected will not remain so. Tracking the evolution of corresponding pairs allows us to 
distinguish different types of magnetic flux evolution. 

Magnetic flux passing through region II exhibits a slippage-like behaviour. Initially connected plasma 
elements above and below the non-ideal region will change their connections as the flow transports the 
magnetic flux linking them into the non-ideal region. On leaving the shadow of the non-ideal region the 
initially connected elements are both transported in the same direction by the flow and a new ideal connection 
is again established for each plasma element. Although this connnection will not be with the initial partner, 
it will be with a plasma element that was initially close to that partner. 

Magnetic flux passing through region III exhibits the type of behaviour most similar to that shown in 
classical 2D reconnection. Again, initially-connected plasma elements above and below the non-ideal region 
loose their connections as the magnetic flux linking them is transported into the non-ideal region. However, 
on leaving the shadow of the non-ideal region the initially-connected plasma elements above and below the 
non-ideal region are transported in different directions by the flow, along opposing 'wings' seen in Figure [H 
and their separation will therefore increase in time, as in the classical 2D reconnection picture. The new 
ideal connection for a plasma element initially above (below) the non-ideal region will be with a plasma 
element that was initially below (above) the non-ideal region in the distant opposing wing. 

We have been able to make a direct comparison between the pure solutions described in Paper I and the 
composite solutions described here since, for our particular choice of vi, the curl of the inertia term in ^ 
vanishes. In principal we could have chosen other ideal flows for this directly comparable analysis that also 
have a curl-free inertial term. One example is that which results from defining the scalar function as the 
function of field-line coordinates given by 
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Figure 5: Separatrices of Wi„ (grey) and Wout black when vi is given by (fT5|) . The same three types of 
reconnective behaviour as Figure [3] are present 



from which we obtain 



-2y'o 
&oA2 



{yx + xz) . 



(13) 



This is also a stagnation flow, but it differs considerably from the flow considered in the previous section; it 
does not cross the separatrices of the projection of Bq onto the a;?;-plane, and is independent of the third 
coordinate, z. When superimposed on the particular solution, however, the same three regions of differing 
flux evolution are present, as illustrated in Figure \E\ The inflow and outflow channels bounded by the 
separatrices of the quasi-flux velocities are now centred around the separatrices of Bq in the z = plane. 
This demonstrates one of the crucial differences between 2D and 2.5D reconnection and the 3D case. The 
crossing of the separatrices by the flow is only a criterion for reconnection in the 2D case. In 3D the difference 
between Wi„ and Wout is the crucial property for reconnection. 

Another example in this class of flows which can be used to form composite solutions is the rotational 
ideal flow arising from the choice 



This ideal flow is rotating in the same sense for all z, and so does not have the effect of bringing flux into 
and away from the non-ideal region. 

For the three flows examined in this section, the reconnection rate, as determined by the integral of 
the parallel electric field along the reconnection line, is identical, but the magnetic flux evolution quite 
different. The distinct types of reconnective behaviour illustrated here, and in Paper I, may be distinguished 




(14) 
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by considering the associated internal and external reconnection rates, as introduced by Hornig (2007). 

jff^ j(f.cxtcrnal j^intcrnal 

dt dt dt 

The external reconnection rate measures the rate at which ux is transported into (and, equivalently, out of) 
the non-ideal region. This rate is always less than or equal to the total reconnection rate. The difference 
between the two is the internal reconnection rate which is non-zero only if there is a circular flow around 
the reconnection line rather than a stagnation flow. 

The stagnation flow examples illustrated in Figures [4] and [5] both correspond to a purely external re- 
connection rate. In this situation the separatrices of the flow (which divide regions II and III) pass though 
the origin, and so the difference in the electric potential between them is equal to the total reconnection 
rate, i.e. the difference in electric potential across the non-ideal region. For an ideal rotational flow, such 
as that arising from the electric potential given by p^ . the external reconnection rate vanishes and the 
reconnection is internal only. Similarly if the ideal flow is zero (as in the case of the particular solutions of 
Paper I) then the reconnection is purely internal. Thus the interpretation of reconnection rate in this way 
allows for a clear distinction between the different types of solutions considered here. 

We note also that a combination of internal and external reconnection is not excluded in these solutions, 
and is expected if a smooth transition between the purely external reconnection solutions illustrated in 
figures m and [5] and the purely internal reconnection found in the particular solution is to be made. Such a 
solution exists when the magnitude of the ideal flow vi is decreased to be the same, or less than, that of the 
non-ideal flow Mpy2- In addition to the three regions of differing space flux evolution described above and 
illustrated in Figure [H the magnetic flux in these mixed solutions would show rotational dynamics within 
part of the HFT. 



4 Accelerating Stagnation Flow 

In a realistic situation we would expect to see a plasma flow that results in a net transfer of magnetic 
energy to kinetic (bulk) energy of the plasma since magnetic energy is the main source of energy in the 
solar corona. This property must be explicitly prescribed here since the model does not include the time- 
dependent processes external (and possibly prior) to the reconnection process that lead to the build-up of 
a current sheet and corresponding plasma flows. Instead these properties are represented in the model via 
boundary conditions on the flow, magnetic field and pressure profiles. Thus in the expansion scheme we may 
ensure an increase in kinetic energy occurs in the reconnection process by requiring vi • j2 x Bp — Vi • V'p2 to 
be, on average, positive over the volume (which is not the case for the above stagnation fiow). This increase 
in kinetic energy may be the result of a transfer of magnetic energy (due to vi ■}2 x Bq), a transfer of thermal 
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energy (due to Vi • V'p2) or a combination of both effects. Numerical experiments suggest considering a 
plasma flow that sharply changes its direction toward the outflow region; in such experiments fast jets of 
plasma emerging from the reconnection region are observed. We examine in this section an ideal plasma 
flow, Vi, which possesses these properties, using a method similar to that of Section (|3.ip . Now, however, 
since the curl of the inertial term in ^ does not vanish, there is a much larger degree of coupling between 
equations ([T]) and ^ , and the particular solution of Paper I can no longer be used as a solution to ^ . 

Just as in the previous section we will consider incompressible solutions; a plasma flow with a faster 
outflow velocity than inflow velocity must have an associated outflow channel that is narrower than its 
inflow channel. To achieve such a flow we impose a non-symmetric function, 0o, for the lowest order electric 
potential, and then deduce the plasma velocity vi from ([T]). For example, we may impose (f>Q as the function 
of field line coordinates given by 

00 = -■^yotanh(xo) , (15) 

and use the inverse field line mappings to deduce an equivalent expression in terms of x, y and z. An 
analytical expression for vi (which is too long to be shown here) is found using ([5]). In the central region 
the flow has a stagnation structure, as shown in Figure El with single inflow and outflow channels that are of 
different widths. Thus, depending on the direction of the flow, and since it is incompressible, an acceleration 
or deceleration of the plasma takes place. The physically relevant case corresponds to the choice ipo > 0, for 
which the outflow direction is the narrower channel along the y-direction, and so the plasma is accelerated 
during the reconnection process. 

Turning now to the lowest-order momentum equation in the expansion scheme, ([3]), we integrate along 
the field lines ([5]), starting from the plane 2; = to deduce the pressure p2' 

P2{x,y,z)^- (vi • V vi) • Bo ds -I-P2 (a;o,yo) ■ (16) 

Js=a 

We first examine solutions obtained when the free function p2 (a;o,yo) is set to zero. Later in the section 
we shall consider another particular example where p2 (xQ^yo) 7^ 0, and show that the choice of this free 
function has a considerable effect on the reconnection process. 

An example of the pressure profile in the case where p2 {xq ,yo) = is shown in Figure[71Ja) . The expression 
obtained for p2 is dependent on ipQ, and so the pressure profile is independent of the flow direction. Thus for 
the case ipo > which we are considering here, a pressure gradient exists along the outflow direction which 
is in the direction of the flow, and so acts to accelerate the plasma. 

The perpendicular component of the current, j2±, can be determined analytically from Q once the 
pressure is given: 

. ^ (-VV2-(vi-V0vi)xBo 
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Figure 6: Stagnation-point structure of the velocity field Vi corresponding to (fiQ — ~ipQyotanh{xo) / A'^ in 
the plane z — 0. As indicated by the closeness of the contours, the plasma has a greater velocity along the 
outflow direction. 




(a) (b) 



Figure 7: (a) The pressure profile p2 and (b) the Lorentz force j2 x Bq in the plane z = for an accelerating 
stagnation fiow vi. The free function p2 {xq, j/o) in equation (|16p has been set to zero. The outflow is aligned 
with the y-axis, and corresponds to the channel of decreasing pressure. 
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A Lorentz force is present within the outflow channels, and is directed away from the central line of minimum 
pressure, as shown in Figure [Tl^b). Thus, since it is not aligned with the flow direction, this force does not 
act to alter the plasma velocity; in this example it is only the pressure gradient which accelerates the 
plasma, causing the fast outflow jets. The quantity vi ■ j2j^ x Bq — vi • V'p2 is, on average, positive over the 
region provided ipo > 0, i.e. provided the flow is accelerated from its inflow to outflow direction. This net 
acceleration, a consequence of the pressure gradient, results in a net transfer of thermal energy to kinetic 
energy. 

The full form of the current j2 may be determined by finding a scalar function A (x, y, z) such that setting 
h = h± + ensures the current is divergence-free. Taking V' • + V'A • Bq = and integrating along 
the field lines gives 

Hx,y,z) = - / V •j2^ (is + A (.To, yo) (17) 

= Hx,y,z) + \{xo,yo), 

where A {xq, yo) is a function that we are free to impose on the solution The current is then given by 

h = {i2± + Hx,y,z)Bo) + X{xo,yo)'Bo = j2 + Js, 

where ^2 is solely determined by the free function A {xq, yo)- The term also determines the current along 
the z-axis because, due to the vanishing divergence of j2j^ along the z-axis, the z-component of j2 vanishes 
there. Equation ^ then implies that the reconnection rate will be determined by this free hmction (together 
with the form of f)), rather than by the ideal flow vi. i.e. governed by the lowest order non-ideal solution. 

As previously mentionned, we can make use of the freedom to choose the function p2 {xq, yo) that arises 
in the integration for p2, given by equation (|16p . The form taken for p2 {xq, yo) will alter the pressure profile, 
and consequently also the current j2. Thus the various choices correspond to additional 3D reconnection 
solutions with differing physical motivations. In particular, a form for p2 {xq, yo) may be imposed such that 
the acceleration of the plasma is driven by the Lorentz force, j2 x Bq, rather than by the pressure gradient, 
— V'p2- One such example is obtained by adding the additional function given (in terms of the fieldline 
coordinates) by 

P2{xo,yo)^P2o e-''^ (yl-xl) (18) 

to the pressure p2- The resultant pressure profile is illustrated in Figure [5]^ a) , where it is seen that the 
pressure gradient is still directed along the outflow channel, but now acts against the direction of the flow. 
Therefore, the acceleration of the plasma is driven entirely by the Lorentz force, with magnetic energy being 
transfered to kinetic energy in the reconnection process. The Lorentz force is illustrated in Figure [8ljb) , 
where the perpendicular component of the current j2 has been deduced using the same method as described 
above. The divergence of j2± along the z-axis remains zero with the inclusion of the extra factor in the 
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(a) (b) 

Figure 8: (a) The pressure profile p2 and (b) the Lorentz force j2 x Bq for the accelerating stagnation flow, 
now with free function p2 (xq, tjq) in equation (|16p given by (|18p . The outflow is aligned with the y-axis, and 
corresponds to the channel of increasing pressure. This figure may be compared with Figure [71 

solution for p2. Thus the reconnection rate in this case is still determined completely by the free function 
A {xotUq), s-nd may therefore be the same as in the previous pressure gradient driven model. 

In both of the examples in this section we must impose a localised resistivity in order to ensure a localised 
non-ideal region. In principle the remaining quantities could then determined numerically using the iterative 
scheme outlined in Paper I. 

5 Conclusions 

In this series of papers we have focussed on the fundamental questions of existence and parametric dependence 
of steady-state 3D MHD magnetic reconnection. Through the use of an expansion scheme, a number of 
examples have been provided that illustrate the very wide range of solutions possible in 3D, a key feature 
of 3D reconnection. At each level of the expansion several free parameters and functions appear, with the 
various possible choices for each of these corresponding to physically different 3D reconnection events. 

One aspect of the diversity in 3D solutions was suggested by the kinematic analysis of Hornig & Priest 
(2003) where it was noted that in 3D, for a given magnetic field, Ohm's law can be decomposed into non-ideal 
and ideal components, the non-ideal solution being termed 'particular' and the sum of the two 'composite'. 
The analysis presented here demonstrates a similar decomposition can be made in a fully dynamic situation 
with the non-ideal and ideal solutions being linked together through the equation of motion. 

In Paper I particular solutions were analysed in a setup with the non-ideal region localised in all 3D and 
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the basic state an X-type current-free equilibrium. The plasma flow in these solutions is purely rotational 
and confined to be within the HFT consisting of all the field-lines threading the non-ideal region. In Section[3] 
of the present work we considered the extent to which the coupling between the non-ideal and ideal flows 
through the equation of motion restricts the form of the ideal solution and examined the effect of the ideal 
solution on the reconnection rate, evolution of flux and energetics. A particular class of ideal flows (for 
which the inertial term can be expressed as a gradient), may be imposed on the particular non- ideal solution 
without altering the form of the current, parallel electric field or (in consequence) the reconnection rate. For 
these flows the coupling between the two solutions is relatively weak, affecting only the pressure term in the 
non-ideal solution. A wide range of flows, both in strength and, more importantly, in profile, belong to this 
class of solution. They may be distinguished by their effect on the evolution of magnetic fiux. 

In general, stagnation flows are expected to be present if classical reconnection is to occur, since they 
allow thin current sheets to be built up and so localised non-ideal regions to become established. A variety 
of symmetric stagnation flows, as considered in Section [31 belong to the class of ideal flows that may be 
used to form composite solutions. These flows bring magnetic flux into the non-ideal region from large 
distances and subsequently remove the flux. Magnetic flux threading particular channels in the centre of the 
region shows similar behaviour to typical 2D reconnection, in the sense that field lines brought in toward 
the non-ideal region reconnect with field lines initially far away, and the separation of initially connected 
plasma elements increases in time after the flux has left the non-ideal region. In the same reconnection event, 
magnetic flux passing through other regions of the domain can be seen to undergo a slippage-like behaviour. 
Although the reconnection rate in the particular and composite reconnection solutions is quantitatively the 
same, its physical interpretation differs. The reconnection, which was completely internal for the particular 
solution is now completely external for these stagnation flows. It was also shown that the 2D criterion for 
reconnection, namely that the plasma flow has to cross the separatrices, does not hold in 3D since we can 
construct solutions which do not have this property but nevertheless have a non- vanishing reconnection rate. 

In these symmetric examples there is no net transfer of magnetic energy to bulk energy of the plasma. 
Non-symmetric stagnation flows, however, such as those considered in Section|3]can convert magnetic energy 
into kinetic energy. These ideal flows show highly curved streamlines, with fast jets of plasma emerging from 
the central region. Although a stronger coupling between the ideal and non-ideal solutions is present in this 
situation, we have shown that, just as in the non- accelerating case, the ideal flow itself does not directly govern 
the reconnection rate. Again the choice of different free functions has been shown to correspond to differing 
physical reconnection scenarios; two particular free functions were used to illustrate how an acceleration of 
the plasma may be driven by a Lorentz force or by a pressure gradient. In general, a combination of both 
effects is also possible. 
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